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In this paper we investigate a solution of the Dirac equation for a spin-| particle in a scalar 
potential well with full spherical symmetry. The energy eigenvalues for the quark particle in s^^^ 
states (with n — ~1) and p^^^ states (with n — 1) are calculated. We also study the continuous 
Dirac wave function for a quark in such a potential, which is not necessarily infinite. Our results, 
at infinite limit, are in good agreement with the MIT bag model. We make some remarks about 
the sharpness value of the wave function on the wall. This model, for finite values of potential, also 
■ - - . could serve as an effective model for the nucleus where U (r) is the effective single particle potential. 
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<U ; I. INTRODUCTION 

Q; 

. In 1964 Gell-Mann and Zweig independently proposed a quark model in which all known hadrons were described 
CnJ ' as bound states of only three fundamental particles P, 0|- Each of these spin-1/2 particles, which Gell-Mann named 

, quarks, has their corresponding antiparticles. However, a precise mechanism of the bound states and quark confine- 

(— I mcnt has still not been established. Gell-Mann's model requires that the valid quark states should be colorless and 

] so a free quark is not allowed. Since a free single quark has not yet been observed, the model assumes that quarks 

• appear as bound states. The baryons are bound states of three quarks (or antiquarks) and the mesons are made up 
1^ ' of one quark and one antiquark. As a model, to describe this behavior, one can think of a hadron as a cavity in which 

the quark wave function is confined, and the cavity is surrounded by the QCD vacuum Therefore in a simple 
model quarks confined in a sphere with radius R and the Dirac equation sould be solved within this sphere. 
. From this point of view, the solution of the Dirac equation in a three-dimensional scalar potential, apart from its 
^ ' interesting theoretical aspects, provide useful tools for studying the properties of elementary particles. For example, 
On the MIT bag model j^-igj and its charily invariant versions, such as the chiral bag model [^-Q and the cloudy bag 

• model |10l4l3| . are some models for describing the physics of the nucleon and other baryons. A bag is a region of 
' space in which quarks and gluons are confined, i.e. they are forced by an external constant pressure B, which can 

be fitted using experimentally determined hadron masses, to move only inside the bag. Historically, Chodos et. al. 
have considered the MIT bag model through the Dirac equation (except the bag pressure B) Q. They solved bag 
equations for the massless Dirac fields in three space dimensions. Their solutions are for the special case of static 
spherical boundary. They also computed charge radius and found it to be 1.0 fm. Degrand et. al., in the other case 
of this model, calculated the masses and the static parameters of the light hadrons In the cloudy bag model a 
baryon is treated as a three-quark bag that is surrounded by a cloud of pions. Thomas et. al. investigated the static 
properties of the nucleon within this model [lo| . They found the bag radius to be about 0.8 fm by a fit to pion-nucleon 
scattering in the (3,3)-resonance region. The chiral bag model for the nucleon is a hybrid of quark and meson degrees 
' of freedom, interpolating the two limits of the skyrme model at i? — and the MIT bag model at i? — > oo 'J„ I 1. 15||. 
In the skyrme model mesons acts as gauge particles so that baryons would interact with each other by the exchange 
of mesons [lB|- Skyrmions are the solutions of the field equations. These solutions are solitons and no longer plain 



X 



wave. One may interpret these skyrmions as coherent states of baryons and excited baryons [17[; however, the physical 
interpretation is still not completely resolved. With this model, it is also possible to calculate nucleon masses and 
other particle properties [3 • Both the MIT bag and skyrme model are useful to calculate masses and other properties 
of hadrons. In Ref. 0, Hosaka and Toki investigated the static properties of the nucleon such as masses and magnetic 
moments as a function of R, in both the original chiral bag model and models with vector mesons. The MIT bag 
model introduces many free parameters for energy corrections that could be helpful in understanding the physical 
processes inside the nucleus. 

It is a curious and complex situation to solve the relativistic quantum mechanics problems in a finite potential 
in comparison with the equivalent problem in nonrelativistic quantum mechanics. In the Dirac equation, the wave 
function is continuous, its first derivative is discontinuous and the second derivative has a very large jump, whereas. 
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in the Schrodinger equation the wave function and its first derivative are continuous, but the second derivative has a 
certain jump related to the potential jump. The solutions of the Dirac equation in a S potential exist in the literature, 
see for example [l9l - [2^ . However, for a finite spherical potential well, as far as we are aware of, there is no solution 
in the literature. 

In this paper, we consider the solution of the Dirac equation in a spherically symmetric scalar potential well, which is 
not necessarily infinite. The origin of this scalar potential could be a strong force that binds quarks together in clusters 
to make more familiar subatomic particles, such as protons and neutrons. It also holds together the atomic nucleus 
and underlies interactions between all particles containing quarks (23l . [23 | . Although this model is not compatible 
with the quark confinement at finite potential, it could serve a dual purpose. First, one could observe the evolution 
of the wave functions as C/q — >■ oo, where one could recover the MIT bag model results. Second, for finite values of Uq 
this could serve as an effective model for the nucleus where U{r) is the effective single particle potential emerging from 
the meson exchange of the nucleons. Here we obtain eigenvalue equations for the energies and numerically calculate 
the energy eigenvalues for the s^^^ and p^^^ states. These states are derived normally from the continuity of the wave 
function and imposing boundary conditions on the cavity ]2S\, with different radii R ~ 0.8, 1, 1.18 fm and the quark 
masses m = 0, 1 fm~^. We then, compare the values of the energy levels in each of specific conditions and also with 
the energy eigenvalue obtained in the previous MIT bag model. The relation between energy eigenvalues and the 
radius of cavity and the mass of the quark is considered. We also obtain the Dirac wave function components for 
a quark particle and depict them in figures. As a result, when the strength of the potential is increased, the wave 
function components on the boundary of cavity fall down and it would have the sharper point, which are in good 
agreement with the MIT bag model . Finally, we obtain the magnitude of sharpness of the relativistic wave function 
component when crossing the wall. 

The paper is organized as follows: Section|TT]is devoted to an introduction of the Dirac equation with a central scalar 
potential. In Sec IIIII we calculate the energy eigenvalues for a Dirac particle in a scalar potential with full spherical 
symmetry for the s^^^ and p-^^^ states. Then we discuss about the sharpness of the wave function components on the 
boundary. Finally in Sec IIVI we summarize our results. 



II. SOLUTIONS OF THE DIRAC EQUATION IN A CENTRAL POTENTIAL (SCALAR COUPLING) 

In the Dirac equation, the scalar potential U{x) and the fourth component of a vector potential, Vo{x) are accompa- 
nied by mass m and energy E, respectively. Although Vq (x) is not a vector potential, since it is the fourth component 
of a four vector, it is called a vector potential. The general form of the equation of motion for a spin-1/2 particle with 
these two potential is (in relativistic units, h = 1 and c = 1) 

[a.p + P{m + Uix)) + Vo{x)]4>ix,t)=i^^^ix,t), (1) 

where p = — iV is the three-dimensional momentum operator. In the above equation a and /3 are the 4x4 Dirac 
matrices which, in the usual representation, are given by 

where / is the 2x2 unit matrix. The subscript i can take the values of 1,2,3, and ai are the 2x2 Pauli matrices. For 
a Dirac particle in a spherically symmetric potential field, the total angular momentum operator J, and the spin-orbit 
matrix operator K = {a.L + 1), commute with the Dirac Hamiltonian. Here L is the orbital angular momentum. 
The complete set of the conservative quantities with their eigenvalues can be written as follows: 

so that, 

[H, K] ^ 0, [H; J] = 0, [J2, J3] = 0, [J, K] = 0, (4) 



and 



± J 



for 



(5) 
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Therefor the quantum number k is a nonzero integer number. Given stationary solutions ip'^,; (x,t) = V'?,. (x) e 
we have 



where and are real square- integrable functions, and yj^' and yj^, can be written in terms of the spherical 

harmonic functions with the relevant Clebsch-Gordan coefficients. Then the two coupled equations for the radial parts 
of the Dirac equation with a given scalar potential turn out to be, 

'^^^^''^ ' ^""-./kM = (m + U{r)-E)g,{r) (7) 



dr 



r 



%^ + ^.9.M = {m + U{r)+E)f,{r). (8) 
dr r 

The solutions of Eqs. ^ and ([U for a scalar potential well, U{r < R) — and U{r > R) = Uq, are the spherical 
Bessel functions and modified spherical Bessel functions of the first kind, for the regions r < R and r > R, respectively. 
Therefore for region 1, r < R with U{r) = 0, and we have 

g^ir) = iVj|^|_i(pr) 

- - ■ p , 



and for region II, r > i? with U{r) — Uq, 



for.<0: ^ ^ "^.T^ti^V ( ^ W 



where p — y/E'^ — m? and q ~ y {m + Uq)'^ — E'^, and N, N' , M, and M' are the normalization factors [1, [25l, [26j. 



III. A DIRAC PARTICLE IN A POTENTIAL WELL WITH FULL SPHERICAL SYMMETRY 

In this section we compute the energy eigenvalues for a Dirac particle in a full spherically symmetric scalar potential. 
Reported eigenvalues, for which the massless spin-1/2 field is confined to an infinite spherical potential well of radius 
R = \ fm, are listed in Table HI These values have been computed through an equation of motion and boundary 
conditions 0j. 



state 


Is, , 


2s, , 




2Pi/, 


E (fm"') 


2.04 


5.40 


3.81 


7.00 



TABLE I: Energy eigenvalues for a massless quark in an infinite potential well with R = 1 fm. 

The solution of the Dirac equation for a particle with mass to, which moves in a spherically symmetric static cavity, 
is physically similar to a scalar potential well with the full spherical symmetry. To confine the wave function of the 
particle in a bag, as the vector potential is zero Vb = 0, the depth of the well should be infinite. This method, which is 
applicable as C/q ^ oo, has been developed by the MIT bag model Q. However, here we consider the scalar potential 
which is not necessarily infinite. 
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FIG. 1: Energy levels for s^. states for a Dirac particle in a spherically symmetric potential well. 



A. The energy Eigenvalues for s^^^ states 



The wave functions for s^^^ states (k = — 1) in a finite scalar potential (with no vector potential, Vq = 0) can be 
written using Eqs. (Pa|) and (jlOal) . For two regions, inside and outside the static spherical cavity, we have 



for r < R 



sin(pr) 



pr 



/_i(r)=-7V 



P_ 

m + E 



ji{pr) = - 



Np 
m + E 



sin(pr) cos(pr) 
{prf pr 



(11a) 



(lib) 



for r > R 



-qr 



qr 



Mq 



m + E + Uo 



ki{qr) 



Mqe' 



m + E + Uo 



1 



1 



qr (qr) 



Using Eqs. (jllal) . (jl2al) and continuity of the g-i{r) at r — R one can find M as follows: 

M = 7V-e«^sin(pi?). 
P 



(12a) 
(12b) 

(13) 



From the normalization condition for wave functions {J^^ [gl{r) + fl{r)] r^dr = 1) and Eqs. ^ and after 
some cumbersome calculations, we get the following expression for N: 



N 



R 



R 



Bin (2pi?) / p2 



2p2 2(m + ^)2 4p3 \{m + Ef 



sin^ {pR) 
p^R{m + Ef 



+ 



sin^(pi?) 
2p2 



1 



1 



q {m + E + Uof 



The continuity of the wave function components aX r — R implies 



-(r<R) 
9 



^^{r>R) 
r=R 9 



r=R 



Now using Eqs. (|lip . and ((T5|) . we find the following expression 

1 



E — m 
c 

E + m 



cot (R\/E^ - mA - 



rn + Uo- E 



1 



R{E + m) \l m + Uo + E R{m + Uo + E) 



= 0. 



(14) 



(15) 



(16) 



This eigenvalue equation gives us the energy of the particle in a scalar potential as a function of R and m for the 
states s. 



"1/2 ■ 
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We have numerically calculated the eigenvalues energy for s^^^ states for different values of R and m and plotted 
them in Fig. [TJ It is understood from numerical solutions [see Fig. [T] and Table |ll] that as ?7o — ^ oo for the case of 
m = and i? = 1 fm the value of the energy ground state becomes 2.0428 and for 2s-^^^ state this value is 5.3960, 
which are in good agreement with previously established results [3|(cf- Table As the potential goes to infinity the 
particle is completely confined inside the bag. We have also depicted the eigenvalues of the s^^ states for different 
values of R and m in Fig. [1] and some of the energy eigenvalues are listed in Table HH In Fig. ^ the wave functions 
for the s^^^ states for a specific case {R — 1 fm, m — 0) are shown. We see that, as the depth of the potential well 
increases, the wave function components on the boundary of cavity become sharper, i.e. if the quarks are turned back 
at the edge of the nucleon by a strong interaction, the wave function will be strongly damped in that region which is 
in line with the MIT bag model. 
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FIG. 2: Wave functions for the Is, , and 5s, , states for R — 1 fm and m — 0. Dot-dashed, dashed and solid lines denote 
Uo = 15 fm~^, Uo = 50 fm~^ and Uo — oo, respectively. We see that as the potential well becomes deeper (bag limit) the wave 
function at the boundary fall down. 



B. Energy eigenvalues for p^^^ states 



We derive the wave functions for , states {k ~ 1) using Eqs. (j9b[) and ()10b|) . as follows: 



for r<R : gi{r) ^ N' ji{pr) ^ N' 
hir) = —^JoKpr) = 



sin(pr) cos(pr) 
[prf pr 
N'p sin {pr) 

m + E pr ' 



(17a) 
(17b) 



for r> R : gi{r) = M'ki{qr) = Af'e"'?'- 



qr [qr) 



Mr) . -M'J"-l±^K4,r) ^ -M' + " 



(18a) 



(18b) 



rn + Uo + E ' ' \ m + Uq + E qr 

As for the s^^^ states we can compute the normalization factor AI' from Eqs. (I17bp and (|18bp . and imply the continuity 
condition of fi{r) at r = R, to have 



M' = -iV'!^i±^^e^«sin(pi?). 

m + E 



(19) 



We have also computed N' from the normalization condition to be 



1 R sm{2pR)^ 1 , i? , sin(2pi?) ^ cos{2pR) 



{m + E)^'-2 



Ap 

sin^ {pR) 1 



2i?p4 2p2 



4pi 
1 



2Rp^ 



1 



{m + Ey 2q 



'1/2 



(20) 
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Now using Eqs. ([15]), ^7} and (HH]), we have 

y^{E-rn) (m + Uo+E) (^1 + 



cot(i?v/£'2 - m2) 



^(i; + to) (m + C/o - £;) 



(21) 



This eigenvalue equation gives the particle energy for a given scalar potential as a function of R and to for p^^^ 
states. We have depicted the eigenvalues for p^^^ states, found numerically from Eq. (PT|). in Fig. [3] At infinite limit 
{Uq oo) we see the energy eigenvalues for Ip^/^ ^-'^'^ ^^'1/2 ^^^^^^ ^'"^ 3.8115 and 7.0020 (for 771 = and R = 1 fm) 
which are in good agreement with previously reported values (see Table H]) Q . Again, in the infinite potential the 
particle is completely confined inside the bag. We have listed some of the energy eigenvalues in Table IIIII We have 
also illustrated the wave functions of the p-^^^ states for R — 1.18 fm and m = 1 fm^^ in Fig. |4l 

Comparing the results of energy eigenvalues for the quark particle in the different static spherical cavity radius i?, 
and different quark masses to, [Figs. [T] and |31 and Tables |lT] and IIII] one can see that an increase in R decreases 
the energy eigenvalues. Also, an increase in to leads to the increase in energy. Therefore, the energy eigenvalues of 
a quark particle confined in a static spherical cavity highly depend on its mass and the radius of the cavity. These 
results are in agreement with the relativistic statistical mechanics. The quark particles inside the bag behave similar 
to a relativistic gas, so that the quarks kinetic pressure is equal to the pressure of the gas 0, H^j- 



E(fin-') 



fi=I fm 





























^■(fin') 

IB 



10 15 20 25 30 35 40 45 50 

f/„(fin-') 




FIG. 3: Energy levels for p^^^ states for a Dirac particle in a spherically symmetric potential well. 



C. Additional remarks about the sharpness of the wave function 



In the case of the scalar potential discussed in previous sections, it is also interesting to investigate the wave 
function continuity at r = R for any k. To study the situation, we use Eqs. ([7]) and ([8|) along with the definitions 
Gk (r) = rg^, (r) and F^, (r) — rf^, (r) for r > R, then we have 



(PG^jr) _ k{k + 1) _ 
dr"^ r'' 



G^{r)~ {m + U{r)y 



dU{r) dG^jr) 
dr dr 



TO 



+ E + U{r) r[m + E + U{r)] 



(22) 



dU{r 



Wherever U{r) has a sharp point, 
whose first derivative is a larger infinity. In this case, to compensate for such a large jump, ^-^^ should have the 



has a certain jump. As U{r) goes to infinity we have the Dirac 5 potential. 



same large jump. Therefore, we can conclude that ' has a jump 



dG^(r) 
dr 

We know that {r) is a continuous wave function. Now integrating Eq. \22\ in the small interval [i? — e, i? + ej, 
and taking the limit e o lead to the zero contribution for the second and third terms. On the other hand, we have 



5{r - R)F{r)dr 




r e [a, b] 
r — a or b 
r i [a, b\ , 



(23) 
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FIG. 4: Wave functions for the p^^^ and 5p^^^ states for R — 1.18 fm and m = 1. Dot-dashed, dashed and solid lines denote 
Uo = 15 fm~^, Uo ~ 50 fm~^ and Uo = co, respectively. We see that as the potential well becomes deeper (bag limit) the wave 
function at the boundary fall down. 



and 



dU (r) 
dr 



UoS [r - R) 



Using Eqs. (f23|) and (|24|) we can compute the remaining terms of integration to get 

R+e\ 



lim 



dG^ir) 



dr 



R[m + E + ^' 



(24) 



(25) 



where G'k{R) denotes the mean value of the first derivative of G'k('') on the wall. This relation gives the sharpness of 
the wave function GK,{r) while crossing the wall. In a similar way, one can obtain the sharpness of F^ir) 



lim 



dr 



R-e ; 



kUq 



R[m-E+!^' 



-F^iR). 



(26) 



IV. CONCLUSIONS 



In this paper we have computed the energy eigenvalues for a Dirac particle in a scalar potential with full spherical 
symmetry. Although in finite potential this model is not compatible with quark confinement, however here we observed 
the evolution of the wave functions as Uq — > oo, where we recovered the MIT bag model results. We have found the 
components of the wave function and depicted them in Figs. [2]and|4l As a conclusion energy eigenvalues of a particle 
for any potential value are positive values. For both cases of s^^^ and p^^^ states in infinite potential the particle is 
completely confined inside the bag. It is pleasing that all of our results are in good agreement with that which exists 
in the literature for the infinite potential limit (the MIT bag model). We also have additional remarks about the 
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continuity and sharpness of the wave function. In the case of the Dirac equation with a finite potential, the wave 
function is continuous, however it has a sharp point at r — R. We have calculated the sharpness of the wave function 
for any k, and we see as the potential goes to the infinity (C/q — )■ oo), the wave function becomes discontinuous. We 
saw the results preserved the relativistic statistical mechanics, so that the quark particle inside the bag behaves as a 
relativistic gas. 
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States 


Potential 


m=0 


m=l fm ^ 


m=0 


m=l fm ^ 


m=0 


m=l fm ^ 




Uo (fm-i) 


/J=0.8 fm 


i?=0.8 fm 


i?=l fm 


R=l fm 


i?=1.18 fm 


i?=1.18 fm 




15 


2.4492 


2.9949 


1.9758 


2.5342 


1.6829 


2.2529 


Is, , 

1/2 


50 


2.5218 


3.0625 


2.0225 


2.5772 


1.7166 


2.2835 




CX) 


2.5535 


3.0932 


2.0428 


2.5966 


1.7312 


2.2973 




15 


6.4655 


6.6608 


5.2178 


5.4262 


4.4451 


4.6663 


2s, , 
1/2 


50 


6.6614 


6.8464 


5.3424 


5.5437 


4.5344 


4.7501 


00 


6.7450 


6.9284 


5.3960 


5.5960 


4.5729 


4.7875 




15 


10.2508 


10.3833 


8.2845 


8.4197 


7.0611 


7.2029 


3^1/2 


50 


10.5885 


10.7020 


8.4922 
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7.3417 


00 
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7.2691 


7.4019 
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14.0606 


11.3090 


11.4184 


9.6498 


9.7590 


4«l/2 


50 


14.4868 


14.5694 


11.6192 


11.7094 


9.8620 


9.9594 




00 


14.6706 


14.7509 


11.7365 


11.8251 


9.9462 


10.0423 




15 






14.2597 


14.3819 


12.2134 


12.3104 


1/2 


50 


18.3745 


18.4400 


14.7382 


14.8094 


12.5096 


12.5863 




00 


18.6098 


18.6725 


14.8878 


14.9572 


12.6168 


12.6921 


y eigenvalues (in fm ^) for the s^^^ states for a Dirac particle 


in a spherically symmet 
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Potential 


m=0 


m=l fm ^ 


m=0 


m=l fm ^ 


m=0 


m=l fm ^ 




Uo (fm-i) 


i?=0.8 fm 


Ji=0.8 fm 


R=l fm 


R=l fm 


Ji=1.18 fm 


i?=1.18 fm 




15 


4.5719 


4.7840 


3.6878 


3.9207 


3.1410 


3.3928 


lPl/2 


50 


4.7056 


4.9098 


3.7738 


4.0007 


3.2030 


3.4498 




00 


4.7644 


4.9669 


3.8115 


4.0370 


3.2301 


3.4756 




15 


8.3835 


8.5171 


6.7686 


6.9107 


5.7671 


5.9189 


2Pl/2 


50 


8.6440 


8.7642 


6.9325 


7.0659 


5.8840 


6.0293 




00 


8.7525 


8.8709 


7.0020 


7.1341 


5.9339 


6.0781 




15 


12.1192 


12.2330 


9.8080 


9.9169 


8.3629 


8.4755 


3Pl/2 


50 


12.5457 


12.6315 


10.0621 


10.1569 


8.5403 


8.6434 




00 


12.7042 


12.7879 


10.1633 


10.2567 


8.6130 


8.7149 




15 






12.8059 


12.9075 


10.9389 


11.0341 


4Pl/2 


50 


16.4353 


16.5026 


13.1823 


13.2561 


11.1888 


11.2689 




00 


16.6445 


16.7093 


13.3156 


13.3878 


11.2844 


11.3632 




15 










13.4807 


13.5764 


"'"'Pi/. 


50 


20.3185 


20.3744 


16.2979 


16.3588 


13.8337 


13.8995 




3C 


20.5799 


2().():i27 




i().r)227 


i:i.9r)2."3 


ii.()l()7 



TABLE III: Energy eigenvalues (in fm '^)for the p^, states for a Dirac particle in a spherically symmetric potential well. 



